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Abstract: Although linear trend removing has often been implemented as a routine 
preprocessing step in resting-state functional magnetic resonance imaging (RS-fMRI) 
data analysis, the spatial distribution of the magnitude of linear trend is still unclear. 
Further, it is interesting whether there will be any difference of the linear trend 
magnitude between different resting-states. For the first aim, we analyzed 5 RS-fMRI 
datasets from 5 different scanners (namely Beijing-Simens-3T, Cambridge-Siemens- 
3T, CCBD-GE750-3T, Milwaukee-GE-3T, and Oulu-GE-1.5T). One-sample t-tests 
on the regression coefficient (i.e., the magnitude of linear trend) were performed for 
each datasets. For the second aim, we used 2 datasets in each of which different states 
were compared, one containing eyes-open resting-state (EO-RS) vs. eyes-closed 
resting-state (EC-RS) and the other containing two steady-state tasks, i.e., real-time 
finger force feedback (RT-FFF) and sham finger force feedback (S-FFF) tasks. Paired 
t-tests were performed between EO-RS and EC-RS, and between RT-FFF and S-FFF. 
One-sample t-tests showed that the spatial pattern of linear trend of RS-fMRI time 
series were quite different between different manufactures. The 3T SIEMENS 
scanners showed positive linear trend in almost all part of the brain, while GE 
scanners showed primarily negative linear trend in most part of the brain. Paired t- 
tests showed some differences between paired conditions; differences between EO-RS 
and EC-RS were mainly in cuneus and eyeballs, and differences between RT-FFF and 
S-FFF were found in the thalamus, anterior cingulate gyrus, and right sensorimotor 


cortex. The current study indicated that, while the manufacturer-dependent linear 


trend of RS-fMRI time series were mostly scanner-related noise, the linear trend may 
also be physiological noise (eyeballs) or even physiologically meaningful, especially 


during steady-state tasks. 
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Introduction 


Since the first study by Biswal and colleagues (Biswal et al., 1995), blood oxygen 
level dependent (BOLD) resting-state functional magnetic resonance imaging (RS- 
fMRI) has been widely utilized to explore the spontaneous or intrinsic brain activity, 
as well as to investigate the pathological progress of brain disorders. There have been 
routine procedures for data preprocessing, e.g., slice timing, head motion correction, 
spatial normalization, spatial smoothing, linear trend removing, and temporal filtering. 
Amongst these procedures, linear trend removing is a well-accepted procedure (Bai et 
al., 2008; Zhang et al., 2008; Fox et al., 2009; Qiu et al., 2011). A common practice to 
remove linear trend is detrending (regression) with general linear model (Bandetini et 
al., 1993; Cox, 1996; Zhang et al., 2008; Fox et al., 2009) during preprocessing in RS- 
fMRI. Although it has been speculated that the linear trend of fMRI signal was a 
system noise arising from scanner instability (Huettel et al., 2004), few studies, 
however, have systematically investigated the spatial distribution of the linear trend of 


RS-fMRI signal in the brain. Further, no study has investigated whether, even a small 


part of, the linear trend of RS-fMRI signal is different between different resting states. 


Such difference might be of physiological importance. 


The first aim of the current study was to systematically investigate the spatial 
distribution of the linear trend of RS-fMRI signal across the brain on data from 
various research sites with different scanners. The second aim was to investigate the 
potential difference of the linear trend between resting-state with eyes closed (EC-RS) 
and resting-state with eyes open (EO-RS). We also compared the linear trend between 
two states of continuous performance, i.e., real-time finger force feedback (RT-FFF) 


and sham finger force feedback (S-FFF) (Dong et al., 2012). 


Materials and methods 


Participants and imaging protocols 


We implemented our analyses on six imaging datasets (See Table 1 for details): 1) 
“Beijing EO/EC” dataset with two resting-state conditions (EC-RS and EO-RS 
without fixation) counterbalanced across participants, 2) “Beijing Feedback” dataset 
with two conditions (RT-FFF and S-FFF, counterbalanced) (Dong et al., 2012) (also 
see below for experiment design), 3) “Cambridge” dataset, 4) “CCBD” dataset, 5) 
“Milwaukee” dataset, and 6) “Oulu” dataset. Datasets 1, 3, 5, 6 were from the “1000 
Functional Connectomes Project” and the “International Neuroimaging Data-sharing 


Initiative” (INDI) (http://fcon 1000.projects.nitrc.org). Dataset 2 was scanned at 


Beijing Normal University Imaging Center. Dataset 4 was acquired from the Center 


for Cognition and Brain Disorders (CCBD) at Hangzhou Normal University. The data 
acquisition of each dataset was approved by the corresponding institutional review 


board. All participants gave written informed consent before scanning. 


The experimental design of “Beijing Feedback” dataset 


The data was from a study designed to investigate the brain mechanism of real-time 
feedback (Dong et al., 2012). During {MRI scanning, the participants were asked to 
continuously maintain a pinch force at about 20 cm H20. In the RT-FFF condition, 
the pinch force was shown to the participant in a real-time way on a monitor. But in 
the sham condition, a pinch force video of another participant was shown to the 
participant. These procedures were implemented by a multiple-channel MRI- 
compatible physiological monitor (model MP150, BIOPAC Systems, Inc., Goleta, 
CA). Detailed experimental procedure was described in the original paper (Dong et al., 


2012). 


Data Preprocessing 


Unless otherwise stated, all preprocessing was performed using the Data Processing 
Assistant for Resting-State fMRI (DPARSF) (Yan and Zang, 2010) 


(http://www.restfmri.net) which is based on Statistical Parametric Mapping (SPM8) 


(http://www.fil.ion.ucl.ac.uk/spm) and Resting-State fMRI Data Analysis Toolkit 


(REST) (Song et al., 2011) (http://www.restfmri.net). The first 10 image volumes 


were discarded for scanner calibration and for participants to get used to the 


circumstance. Then slice acquisition dependent time shifts were corrected per volume. 
The time series of images for each participant were realigned using a six-parameter 
(rigid body) linear transformation with a two-pass procedure (registered to the first 
image and then registered to the mean of the images after the first realignment). 
Spatial smoothing was performed using a 4 mm Gaussian kernel. All volumes were 
then normalized with EPI template to obtain transformations from individual native 


space to MNI space. 


Linear trend estimation 


Linear regression analysis was performed in a voxel-by-voxel way to calculate the 
regression coefficient (beta) of the time series against time at individual native space. 
The beta maps were then registered into MNI space with re-sampled 3 mm? cubic 
voxels by using transformations acquired from spatial normalization procedure. One- 
sample t-tests were performed for each dataset on the beta maps. And paired t-tests 
were used to compare differences of beta between EO-RS and EC-RS as well as 
between RT-FFF and S-FFF for “Beijing EO/EC” dataset and “Beijing Feedback” 
dataset, respectively. All statistical t maps were corrected by Gaussian random field 
(GRF, voxel p < 0.001, cluster p < 0.05) within the whole bounding box (271633 
voxels in all) in order to investigate the linear trend both inside and outside of the 


brain. 


Results 


The spatial distribution of the linear trend of {MRI signal 


The linear trend maps were shown in figure 1. The pattern of linear trend of fMRI 
signal showed a manufacturer-specific spatial distribution: SIMENS scanners showed 
positive linear drift, and GE scanners showed negative linear drift in most parts of the 
brain. In addition, linear trend in the white matter is stronger than that in the gray 
matter for all datasets. Specifically, the Cambridge SIEMENS scanner showed 
apparently higher drift than Beijing SIEMENS scanner, partially due to its larger 
sample size. For the 3 GE scanners, the drift was apparently the highest for Oulu, the 
lowest for CCBD, and moderate for Milwaukee. Although the GE scanners showed 
primarily decreasing drift, the CCBD and Oulu scanners showed an increasing drift in 


the orbitofrontal midline area. 


Difference of linear trend between different states 


In the EO-RS vs. EC-RS map, visual area (V1) showed a significantly smaller 
linear trend in the EO-RS state than EC-RS state, and eyeball showed a significantly 
larger linear trend in the EO-RS state than EC-RS state (Figure 2, Table 2). In the RT- 
FFF vs. S-FFF map, the bilateral thalamus, bilateral anterior cingulate gyrus (BA 32), 
and the right sensorimotor cortex showed significantly smaller linear trend in the RT- 


FFF state than S-FFF state (Figure 2, Table 2). 


Discussion and Conclusion 


This study demonstrated a manufacturer-specific pattern of linear trend of RS-f{MRI 
BOLD signal, suggesting a thermal noise over time. However, it is difficult to 
interpret the opposite direction of linear drift between SIEMENS and GE scanners. 
Our results would be helpful for the manufacturers to further analyze the linear drift 


and hence to improve the signal stability. 


Although the above results strongly support that removing linear trend should be a 
routine procedure during preprocessing of RS-fMRI, the significant differences of 
linear trend between EO-RS and EC-RS as well as between RT-FFF and S-FFF 
indicated that the linear trend of RS-fMRI BOLD signal is not completely scanner 
noise. It would be informative to compare the linear trend between groups or between 


conditions before removing the linear trend. 
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Table 1. Detailed information of datasets. 


Dataset N Sex Handness Age EC-RS vs. EO-RS Scanner Magnet TR TE Filp angle Slice Time 
(M/F) (R/L) mean/SD (ms) (ms) (°) number points 
Beijing EO/EC 42” 21/21 42/0 22.3/1.6 Two counterbalanced SIEMENS 3T 2000 30 90 33 240 
sessions TRIO 
Beijing 38° 19/19 38/0 22.3/1.6 Two counterbalanced SIEMENS 3T 2000 30 90 33 240 
Feedback sessions TRIO 
Cambridge 198 75/123 171/27 21.0/2.3 Open SIEMENS 3T 3000 30 85 47 119° 
TRIO 
CCBD 41 21/20 41/0 23.9/1.6 Close GE 750 3T 2000 30 90 43 240 
Milwaukee 45 15/30 * 53.4/5.6 t GE 3T 2000 £ * 64 175° 
Oulu 95° 31/64 86/9 21.6/0.6 Open, fixation GE 1.5T 1800 zi zi 28 245° 


° There were 48 participants in “Beijing EO/EC” dataset and 43 in “Beijing Feedback” dataset. After excluding four participants from “Beijing EO/EC” dataset (due to bad spatial normalization 


or left handedness) and five participants from “Beijing Feedback” dataset (due to technical problems or excessive head motion) (Dong et al., 2012), 42 and 38 participants were selected with 


session order, sex and age for counterbalancing between EO-RS and EC-RS or between RT-FFF or S-FFF, respectively. 


b Accordingly the first 5 time points of each time series were discarded in preprocessing (for more details: http://www.nitrc.org/docman/view.php/296/716/fcon 1000 Preprocessing.pdf), 5 


extra time points were removed in our studv. 


“There were 103 participants in “Oulu” dataset, eight participants were removed because of bad spatial normalization or excessive head motion. 


Table 2. Detailed information between different states in beta maps. 


States Regions BA Cluster volume (mm?) t-score of peak voxel Coordinates of peak voxel 

EO-RS vs. EC-RS Right eye ball N/A 4806 5.3036 42, 63, -39 
Left eye ball N/A 5076 5.6542 -36, 51, -51 
Sub-Gyral N/A 19737 -6.7337 27, 33, 15 
Cuneus 19 810 -4.7407 3, -81, 30 
Sub-Gyral N/A 1161 -4.6723 -27, -66, -9 
Cuneus, N/A 8910 -6.5379 30, -57, -3 
Sub-Gyral N/A 3969 -6.5393 48, -30, -12 
Sub-Gyral N/A 1593 -4.6042 33, -66, 21 
Sub-Gyral N/A 3159 -4.9170 -15, 12,27 

RT-FFF vs. S-FFF Thalamus, Extra-Nuclear N/A 3645 -4.7279 6, -6, -3 


Cingulate Gyrus 32 2457 -4.7617 -6, 15, 36 


Extra-Nuclear 


Extra-Nuclear 


Extra-Nuclear 


Medial Frontal Gyrus, SMA 


Inferior Parietal Lobule 


Precentral Gyrus 


N/A 


N/A 


N/A 


N/A 


N/A 


N/A 
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702 


594 


756 


1755 


621 


-4.6684 


-4.5915 


-4.4860 


-4.7441 


-4.5269 


-5.5698 


-30, 15, 0 


30, -36, 21 


-24, -27, 18 


12, -3, 54 


36, -39, 54 


33, -24, 57 
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Figure 1. Result of one-sample t-tests (p<0.001, cluster's p<0.05, GRF corrected) on the beta 
value acquired with SIEMENS (EC-RS, EO-RS, RT-FFE S-FFF and Cambridge) and GE 


(CCBD, Milwaukee and Oulu) scanners. N is the number of subjects for each dataset. 


Different States 


EO-RS vs. EC-RS 


RT-FFF vs. S-FFF 


Figure 2. Results of paired t-tests (p<0.001, cluster’s p>0.05, GRF corrected) between EO-RS 


and EC-RS (upper row), and between RT-FFF and S-FFF (lower row). N is the number of 


subjects. 


